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A Bethe Ansatz solution of a (modified) Hubbard chain with a Kondo impurity of arbitrary spin 
S at a highly symmetric line of parameter space is proposed and explored. Our results confirm 
the existence of a strong-coupling (line of) fixed-point (s) with ferromagnetic Kondo coupling as 
first hypothetized by Furusaki and Nagaosa on the basis of perturbative renormalization group 
calculations. For on-site Hubbard repulsion and ferromagnetic Kondo exchange, the ground state 
has spin S — 1/2, i.e. , is a singlet when S = 1/2. The contributions of the impurity to the magnetic 
■ susceptibility and low-temperature specific heat are discussed. While the Wilson ratio is unity in 

| the half-filled band, it is found to be a function of density and interaction away from half-filling. 
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The recent developments in the field of high-temperature superconductors and heavy-fermion materials have stim- 
ulated a renewal of interest into the theory of impurities in quantum liquids and of strongly correlated Fermi-systems. 
^ , The search for non-Fermi-liquid fixed points has gavitated around two prototypal models and their respective descen- 
dants: the Hubbard- (t — J, etc.) and the Anderson (single impurity and periodic Kondo, etc.) models. In the last 
decades, these Hamiltonians and their associated families have become the most extensively studied many-fermion 
problems in condensed matter theory. The progress has been tremendous and the literature in the field is considerable. 
For the sake of brevity , we refer to a few reviews pj and the numerous references therein. 

In the present letter, we want to discuss a Hubbard chain with a Kondo impurity that is integrable on a particular 
line of parameter space. Our model emerges as a hybride of the one-dimensional Hubbard model 0] and the standard 
Kondo model §{§: 
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where c\ a creates an electron of spin a =f, j in a Wannier state centered at site i, rii a — c' ia .Ci t(7 ; t denotes the 

t-H ' hopping integral, U the screened on-site Coulomb repulsion, J the exchange coupling the localized impurity spin 

S to the conduction electrons and V a local charge interaction. In equation (|l]) we have included the chemical 

potential fj,, the Zeeman-term with h — g e /iBH and 7 = g s /g e the gyromagnetic ratio for the ion impurity in units 

of the conduction-band electron g e factor. We assume a chain of L sites and impose periodic boundary conditions. 

The model of equation (jl]) ignores important physical ingredients [|| such as orbital degeneracy, spin orbit coupling, 

, crystal-field effects, and so on. 

To our knowledge, the Hubbard-Kondo (H-K) model of equation (|l]) has not been solved previously, though various 

solutions of related models on the continuum and lattice have been proposed (see [j6|dTo|). As is well known the Kondo 

model can be derived from the Anderson model (see Jll[]) as a low-energy effective hamiltonian via a Schrieffer- Wolff 

transformation in the regime where the Anderson model has a local moment and antiferromagnetic exchange with 

the conduction electrons (a potential scattering also arises in the reduction process). We note that for the H-K model 

considered in this work we assume a priori the magnitude of the exchange J to be arbitrary. Furthermore, due to 

the presence of the electron-electron interactions, the H-K model is a genuine one-dimensional model in contrast to 

the original three-dimensional Kondo model, which is solved exactly by keeping only s-wave states, i.e., left-going 

electrons scattered forward by the Kondo exchange operator. The Hamiltonian (Q) inherits various symmetries of 

its progenitors (Hubbard and Kondo) provided appropriate rotations in spin space are performed []l2|| . The most 

remarkable property of the model (when appropriately modified, see below) is however its integrability at J = —U/2 

and V = or V = — ^(S + 1), as we now discuss. 

To solve the eigenvalue problem associated with the Hamiltonian of equation (Q) , we make the following Ansatz for 
the N-fermion problem: 

^...^ (xi, . . . , x N ) = ®o(Xq)Q{xq){-1) P A <jq1 ..., tqn (Q\PQ) exp (i ^ k Pj x^ (2) 
Q,p 
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where Q(xq) = 0(xqi < ... < xqn) — 1 when the ordering associated to the permutation Q of the coordinates is 
respected and otherwise zero; &o(x~) — Y[j ®i v j x j), where Vj{Q) = ±1 according to whether the coordinate is to 

the left or to the right of the impurity as specified by the permutation Q of (N + 1) objects, and where <d(x) = 1 if 
x > 0, and otherwise zero with the prescription 9(0) = 1/2. 

A straightforward substitution of the Ansatz (|J) into the Schrodinger equation H^/ = leads to difficulties which 
however are only apparent. Indeed, we need a prescription of how to treat the electrons as they approach the 
impurity. For the sake of simplicity consider two particles away from the impurity with amplitude 1 i , <Jl ,a 2 ( x i: x 2)- The 
Hamiltonian is modified near the impurity as follows: when, particle one, say, hops from the right onto the impurity, 
the amplitude becomes ^/ CTli(T2 (0 + , x 2 ), while if it hops from the left we obtain '^^(O", x%). Here 0* means that we 
have taken the limit to zero from the right (left) when approaching the impurity from the right (left). The amplitude 
for particle one to be at site x\ — (impurity site) is defined by ^a 1 ,a 2 (0, £2) = [&a 1 ,cr 2 (0 + ,£2) + ^a 1 ,a 2 {0~ , x 2)]/2, 
and so there is no ambiguity as to a particle sitting on the impurity site. Similarly for particle two. Furthermore, due 
to the presence of the impurity, the four amplitudes \E' (TljCT2 (0 ± , ± ) cannot coincide for our ansatz to be consistent. 
Yet, we can define the amplitude for two particles of opposite spin to be at the impurity site as the average of these four 
amplitudes. With these definitions, the eigenvalue problem can be dealt with in the usual way near the impurity while 
at the impurity site additional terms occur, that cancel exactly those coming from the electron-impurity interaction. 
The idea is readily extended to the iV-particle case . That the low-energy long- wave length physics of the original 
Hubbard model with Kondo impurity should not be affected by our strategy can be seen as follows: the model we 
discuss here can be viewed as a lattice regularization of a continuum Hamiltonian solved in Ref. H (see comments 
in the bibliography). We will provide a complete discussion of this point together with numerical computations 
elsewhere Q. 

The Schrodinger equation imposes linear relations between the coefficients A(Q\PQ) in the Ansatz (||) and determines 
the many-body scattering matrix. The electron-electron scattering matrix is identical to that found in the Hubbard 
model (see 0) while the electron-impurity scattering matrix reads 



_ , . . , 2itsmkj - (V + JPi s ) 

ILJsmfc,) — 

2itsmkj + (V + JP^ 3 ) 
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where s'mkj is the rapidity of the colliding electron, Pi s = | + er ■ S and V — V — J/2. The many-body S-matrix 
factorizes into two-body scattering matrices provided the consistency conditions (Yang-Baxter equations) are fulfilled. 
Since the latter are identical in form to those occuring in the Kondo model, we refer the reader to the literature (see 
for example f|,fl|]). The Yang-Baxter equations require then J = — ^ and V = ±¥(£ + 1/2). 

We can express the periodic boundary conditions in terms of the two-body scattering matrices, TjA ai ... aN s{I\I) = 
e iL A ai ... aN s(I\I), where / denotes the identity permutation and Tj = Sjj+\ . . . SjNRjsSji . . . Sjj-i represents the 
transfer matrix that carries a particle around the ring. The eigenvalue problem can now be solved by the technique 
of the monodromy matrix (Quantum Scattering Method, see for example @,||). For V = ±f (S + 1/2), we find the 
following set of coupled algebraic equations : 



z 2 



_ -py Afj - A a + ig _ A a - igS yr Oj - A a + zf 
A p -A a - ig A Q + igS AA ay - A a - if 

where N c = TV is the number of charges and N s the number of down spins in the system, ou — smkj (j = 1. . . N c ), 
denote the charge rapidities, A Q (a = 1 . . . N s ) the spin rapidities and g — U/2t. Equations (0) resemble those derived 
by Lieb and Wu for the Hubbard model, yet include an additional phase shift due to the presence of the impurity. 
Taking the logarithm of equations (|J) leads to a system of coupled transcendental equations parametrized by two sets 
of quantum numbers I" (a = c, s), which are intergers or half-odd integers depending on the values of N Cl N s and 
L. Notice that the impurity shifts the quantum numbers of the Hubbard model and that in contrast to the Kondo 
model, the electron-impurity phase-shift depends on the state of motion of the charges and spin degrees of freedom 
through their rapidities. In the thermodynamic limit, we introduce the distributions p c {k) and /5 S (A) of the "charge" 
and "spin" rapidities [0] that obey coupled integral equations of Fredholm type. Since the method is by now standard, 
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we refer the reader for the details to [12] and discuss here the physics of the impurity. 

In this letter we focus on the contribution of the impurity to the magnetization and the low-temperature thermody- 
namic properties of the system. The magnetization in zero-field h = is given by M — Alh + j;Mj with Mh = and 
Mi = 5 — 1/2, i.e., the ground-state is 25-fold degenerate and the spin of the impurity is partially screened. When 
5 = 1/2 the ground-state is a singlet despite the ferromagnetic Kondo exchange coupling. This result is somewhat 
counter-intuitive, especially in view of a theorem due to Mattis [S3 which states that for a ferromagnetic Kondo model 
in the absence of electron-electron interactions the ground-state has spin 5 + 1/2, i.e., is a triplet for spin 5 = 1/2. 
Our results can however be understood if we remember that as soon as U is finite, we scale to a (line of) strong- 
coupling fixed point (s): the screening of the local spin is a consequence of the presence of relevant (with respect to 
the Luttingcr liquid fixed point) local backscattering at the impurity site as seen from bosonization of the H-K model 
jL5[|L6). The existence of a strong-coupling fixed-point has tentatively been suggested by Furusaki and Nagaosa pl| 
using perturbative renormalization group techniques. The ferromagnetic Kondo exchange coupling makes the electron 
spin on the impurity ion to align with the localized magnetic moment while the neighboring electrons rearrange in 
such a way as to screen the resulting enhanced local magnetization. 

At arbitrary magnetic fields, the problem has to be solved numerically |l2|. In the following, we consider the weak 
magnetic field regime which can be treated analytically. Let us first discuss the half- filled band n c — 1. Because 
umklapp-scattering is relevant at half-filling, this case was not discussed in ref. Jl^|. In the presence of a weak h > 0, 
the bulk magnetization is (see for example p8|jlS| ]) M/, = j° , where I u (z), for v = 0, 1, denote the Bessel func- 
tions of the first kind of imaginary argument, and the contribution of the impurity is obtained by the Wiener-Hopf 
method 11711 as 



M t « 5 - 1 + -ij J™ ^ sin [(25 - 1)ttw] T Q - cu^j exp (2w Hh/T H )) 

where g\ \sx{h/Tjj)\ ^> 1 and where we have set 7=1. Note that formula (|5|) resembles the one found in the 
antiferromagnetic Kondo model. Moreover, in equation (||) the coupling constant g enters through the expression 
\a.{h/Tji)- The magnetic scale Th (in units of ks Bolzmann) has been defined by Th = (2-7T / 'e) 1 / 2 ?^, where "the 
Kondo-Hubbard temperature", Tk = 4tli(n/g), represents a strong coupling scale controlled by the band-width of 
the conduction-band and a Bessel function that interpolates betwen the strong- {g — ► oo) and weak- (g — * 0) coupling 
limits of the Hubbard model. We expect that for T -C Tk all physical quantities depend only on the energy scale Tk- 
The case (U = —2 J < 0) of attractive on-site interaction and antiferromagnetic exchange is also interesting and will 
be discussed in fl^| . 

We distinguish the compensated (5 = 1/2) and uncompensated (5 > 1/2) cases. When 5 = 1/2, the magnetization 
can be written as 



i^- (-D' n + i/2\ a *" 2 > i h \ (2M > 

^^w^{.) he) (6) 

where h -C Th- To lowest order in the field the magnetization is linear in h, i.e., the impurity magnetic susceptibility 
is a constant Xi ~ ^-I^Tk- This is consistent with the picture that the impurity spin is screened at zero temperature 
by the conduction spin carriers. Here, Tk represents the "binding energy" of a complex of electrons and the impurity 
spin. The strong-coupling limit of the Hubbard model g 3> 1 yields an enhanced susceptibility Xi ~ 9 727r 2 £ identical 
to the bulk susceptibility |is| per site \h ~ U/4ir 2 t 2 , i.e., at weak but finite field the electrons are less and less effective 
in screening the impurity spin with increasing g — U/2t = — J jit — > oo. The zero-temperature susceptibility diverges 
when g — > oo as expected for decoupled spins. Let us recall that in this limit the model maps (via second order virtual 
processes) to a Heisenberg chain with an antiferromagnetic effective exchange J e ~ At 2 /U and an electron spin trapped 
at the impurity site. A simple picture is that of a complex of four spins built from a doublet of the impurity spin 
and two neighboring electron spins coupled via antiferromagnetic effective interaction to an electron spin sitting on 
the impurity site. In the limit of weak coupling, jCl, the magnetic susceptibility behaves as Xi ~ iuflg ex P \ ~a) 
while the bulk susceptibility \h ~ (27ri) _1 is finite, i.e., the impurity susceptibility is strongly suppressed. Here, the 
naive (mean-field) picture of a spin-density wave coupled to a local magnetization applies: the former has effectively 
swallowed the impurity spin. While in a singlet ground-state, we can pass continuously, by varying the Hubbard 
repulsion (Kondo exchange) U = —2J, from a regime where the impurity spin has disappeared locally to a regime 
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where the impurity site has decoupled from the remaining spins. 

For S > 1/2 the dominant terms in the impurity magnetization reflect the physics of a strong coupling regime, 



Al, 



(5-1/2) 



1 \n(\n(T H /h)) 
HT H /h) 2 [\n{T H /h)f 



(7) 



The corrections to the above formula are exponentially small. The resulting susceptibility is singular in the limit of 
small field, a result that can be interpreted in terms of a reduced spin (S — 1/2) having a ferromagnetic interaction 
with the conduction-band electrons. 

It can be shown jl2| that at half-filling (as in the standard antiferromagnetic Kondo model) there is a duality on 
substitution of S by (S — 1/2) between the low (T <C Tk) and high (T ^> Tk) temperatures as seen for example in 
the magnetic susceptibility 



Xi(T) 
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The low temperature specific heat at half-filling for spin S — 1/2 has contributions only from the spin degrees of 
freedom, i.e., C v — C% + C\/L — [l + p l s (q s )/ Lpg(q s )] , where the spinon velocity v s = 2tI\(-K / g) / Iq(-k / g) and the 
bulk specific heat = Iq(it / 'g) /2wtli(ir / 'g) (see see for example J20) ) . For S = 1/2, we have C % v /C^ = \/I {-n/g), 
i.e., Cl/Cy(g — » oo) w 1 and C l v /Cy(g — ► 0) « n^/2~g exp ^— ~Y again a manifestation of the fact that the itinerant 
electrons become more effective in screening the impurity in the limit g — > 0. At half- filling and for spin S = 1/2, the 
Wilson ratio is therefore : B/ 1 / 2 ^ 1 — *\ ^ h h = 1. Near the saturation field h c (see for example the susceptibility 

and specific heat behave as Xi/Xh = CJ/Cj ~ ^/SgF(g) where F{g) is a function of g. We thus infer R( s>1 / 2 ') = 1, 
for h near h c , as in zero-field. In the weak coupling limit, we have Xi/Xh = Cl/Cy « 1/gS while in the strong 
coupling limit, the above ratio becomes independent of the coupling constant g , i.e., Xi/Xh = C\/C } ^ « 1/25. 
We now consider the strong coupling limit at arbitrary band filling. For a compensated impurity, the magnetization 
at weak fields is Mj ~ 2h/?:h c , which is proportional to the bulk magnetization so that Xi/Xh ~ l/n c . The specific 
heat ratio simplifies in this limit and the Wilson ratio becomes 



1 + 



2g sin(7rn c ) 
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ergo is not a universal number, i.e., is a function of the interaction g and density. The above expression is valid 

provided 8 = 1 — n c 3> l/.g. 

At small density n c <C 1 we can estimate the specific heat and susceptibility ratios for S = 1/2. We find a Wilson 
ratio R^ 1 / 2 } w 1 + K 2 n c /g, again a non-universal number as in the large g limit. For spin S > 1/2 and away from 
half- filling, the magnetization is identical in form to that found at half- filling on the substitution Tu — > TS, where 
the new energy scale reads « (2-k / e) 1 / 2 (2t/ g)(Tm c ) 3 . The physical interpretation of this result is similar to that 
of the above case, except that now the band width of the conduction electrons scales to zero with the charge density. 
In conclusion, we have solved and discussed a Hubbard chain with repulsive electron-electron interactions and ferro- 
magnetic impurity exchange. Our solution confirms the existence of a (line of) strong-coupling fixed point (s), which 
had been tentatively suggested by Furusaki and Nagaosa on the basis of perturbative renormalization. 
The authors are grateful to J. -Ph. Ansermet, C. Gruber, P. Nozieres and T.M. Rice for discussions. 
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